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1. 0 ABSTRACT 



nMRSIM(E) ' is designed to calculate theoretical spectra of compounds contain- 
ing nuclei of spin one-half , including hydrogen, fluorine, and carb6n-13. 
Chemical shifts for each nucleus and coupling constants between nuclei are 
input from the Teletype* Calculated line spectra are displayed on the VR12 
scope. T^/zenty-f ive calibration points are displayed across the X axis. Out- 
put is to LINCtape and high- or lovz-speed paper tape punch, as well as to 
the TelotypG. NMRSIM has options that allow varying all the parameters as 
well as offsets for enhancement of resolution. Chemical shifts and coupling 
constants may be adjusted continuously until the displayed theoretical spec- 
trum is acceptable. Spectra may be read back and displayed from LINCtape or 
paper tape; several psectra may be meryed in this mode of operation, thereby 
allowing the simulation of large spin systems or mixtures of compounds. 

2.0 EQUIPMENT 

A PDP-12A equipped with KW12A real time clock and MC12 4K Memory Extension is re- 
quired. 

3.0 OPERATING PROCEDURES 

NMRSIM (E) is loaded from a LAP6-DIAL^ tape containing the program by the command 

LO NMRSIM (E) ,0 ) 

Refer to the LAP6-DIAL Programmer's Manual, DEC-12-SE2B-D for further details. 

The restarting procedure is as follows: 

1) Press the STOP key on the PDP-12A. 

2) Press I/O PRESET. 

3) Press START. 

NMRSIM then returns to the display mode. 
1 

NMRSIME utilizes the EAE option for the PDP-12. 

2 

LAP6-DIAL is hereafter referred to as DIAL. 



4.0 SPECTRUM GENERATION 



After the program is loaded, it prints a series of messages on the Teletype 
to specialize the parameters for the experiment. The user types a .reply to 
each message and terminates the response by pressing the RETURN ke^-. If an 
illegal response is typed, a ? is printed on the Teletype and the rr.essage'is 
repeated. The RUBOUT key can be used to erase incorrectly typed characters 
before a terminator is typed. The messages are listed below with their ac- 
ceptable responses. Refer to section 8.0 for a sample dialogue. 

4 . 1 COMMENTS : 



Any amount of commentary can be typed at this time using any of the Teletype 
keys. After entering appropriate titles, type CTRL/A to advance to the next 
message. 



4 . 2 WANT PAPER TAPE I/O? 

Only a response of Y or N is acceptable. A reply of N stores the position 
and intensity of all calculated transitions on LINCtape only. A reply of 
Y outputs the position and intensity of all calculated transitions on paper 
tape in addition to LINCtape. After a reply of Y, a second question is 

printed : 



WANT HIGH SPEED READER-PUNCH? 
Type Y to use the high-speed punch, type N to output to the Teletype punch. 
4.3 NO. OF SPINS = 



Enter the nximber of spins; values from 1 through 6 are acceptable. 

4.4 OFFSET & WIDTH 

These two parameters define the range of the X axis on the scope. OFFSET 
is the value in Hz of the first point rinplayed along the horizontal axis. 

WIDTH is thr, total range in llz of t.ho horizontal axis. It is important 
to kr;cr^ in mind that the axis will bod splayed from right to leFt corres- 
ponding to th(.t standard display mocJo for NMR dai.ci. 
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4 . 5 CHEMICAL SHIFTS : 



Enter the value m Hz of the CHEMICAL SHIFTS in order. Type a comma to 
separate values: 

1, 2, . , . , n 

4 . 6 COUPLING CONSTANTS 

Enter the value in Hz of the coupling constants in order: 

2, ^1 3, ^^1 4, ^In, 3, ^2 4, ^2n, ^Cn-l)n 

Sjjecify every coupling constant required, even if the value is zero. 
Type a comma after each value. 

4 . 7 BLOCK, U: 

Type the starting block in octal notation, a comma, and the number of 

the LINCtape transport unit to be used to store the peak parameters of 
the calculated transitions. The number of blocks required is dependent 
upon the size of the spin system desired. 

The computer now calculates the theoretical spectra and automatically 
stores the data on the LINCtape. 

5 . 0 COMMANDS 

When dis[. laying a line spectrum the Teletype is active and can be used 

to issue comrnrinds defining the next task for the computer. Six options 

are available to the user at this time. Each of these commands is issued 
by typing the appropriate letter follov/ed by a colon. The command may be 

changed by using the RUBOUT key before the colon has been typed. The displ 
can be scaled v/ith Right Switches 9, 10 and 11. 

3. 1 R (RESTART) 

Restart prints the message COMMENTS (section 4.1) and restarts the pro- 
gram. It is used to calculate a complete set of new data. Note that 
paper tape status can not be changed after the initial responses to the 

messages . 



5.2 C (COUPLING CO'I^STAMTS ) 



Returns to section 4,6 and requests nev/ parameters. Used for open loop 
iteration to find the best coupling constants. 

5.3 O (OFFSET h< WIDTH) 

This command allows a new section of the X axis to be displayed on the 
scope so that any section of the calculated spectrum can be displayed 
on the scope. The OFFSET option is useful for resolution enhancement. 

To exercise this. option, there must be NMR data on LINCtape or paper 
tape (prepared by the P: command). If paper tape was requested during 
initialization (section 4.2), the question 

TAPE INPUT? 

is asked. Type Y to road data from the high speed or Teletype reader and 
v/rite it on LlllCtapc before display. Type N to read data from LINCtape. 

^ , 4 D (DIAL; 

Program returns control to the DIAL Operating System. 
5 . 5 L (LIST J 

The L command lists intensity and energy data on the Teletype in tabular 
format. After typing L, the message 

MIN. INTENSITY: 

is printed. Type an integer number betv^een 0 and 1000. If the intensity 
of the transition is greater than the entered value ^ the transition will 
be in the output list. This question thus serves as a threshold on the 
intensities that are listed, A second question 

. DISPLAY LIST 

is thr;n priritrjr]. Typo Y Lo output (.)ri ly the lines wh. ich fire di.s[)layed on 

the scope. 'fype U to output all linos in the tta}^le» There are usually 

many more lines than are displayed on the scope because the scope cannot 
resolve closely adjacent lines. 



1 

If the Teletype punch is turned on during a List operation ^ a binary 
coded decimal paper tape is produced v;hich may be used as inpu.t to CATACAL 
(DEC~12--UW1A-D) to produce a spectrum consisting of peaks instead of lines. 
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5.6 P (PUNCH) 

The accumulated spectrum can be punched on the high-speed reader/punch 
by issuing the punch conimand. 

6. 0 OPEIWriMG SUGGESTIONS 

A useful starting point for many 60 MHz NMR (hydrogen) spectra of organic 
molecules is a sweep offset of 0 Hz and a sweep width of 500 Hz, For best 
results, however, all peaks should eventually be examined at a considerably 
smaller sv/eep width, since rounding errors may give misleading scope dis- 
plays with large sweep widths* It should be noted that increasing the sweep 
offset moves peaks to the right, and decreasing the sweep offset moves peaks 
to the left. Once the peaks are observed, the proper sweep offset and width 
can be determined from the 25 calibration marks placed across the X axis. 

For purposes of estimating a logical minimum intensity, remember that the 
calibration on the scope on the Y axis is always 0 to 1000. 

1 . 0 NMR BACKGROUND 

The discussion below assumes that the nuclei being studied are hydrogen, 
fluorine, carbon- 13 or other nuclei of: spin 1/2. 

Each nucleus in a rriolocule can be assigned a spin of ei.ther a or 3 . For 
a syHt^:m containiny n nuclcj i. , there are spin states |)Ossible. Each of 
these states is called a basic product function. The energy of a basic 
product function can be calculated from equation 1 belov/ where s^ = -1/2 
for a spin and + 1/2 for a spin, is the frequency of absorption in 

Hz, Tj^j- + 1/4 if i and j have the same spin and - 1/4 if i and j have op- 
posite spin, and J\. ^ is the coupling constant between i and j in cps. 

n 

1 H = (s . u. + >" , T. . J. , ) 

uu 1 1 3>i 13 

The energy of interaction between two basic product functions is 0 unless 

both have the same number of a and B spins « V^hen this condition is ful- 
filled, the interaction energy i.s calculated from equation 2 where U=l if 
the basic product functions differ only in the interchange of spins i and 
j (J.j as de-fined above). Otherwise.!, \y-^0. 

2 U - 1/2U] . . 

uv i. 3 
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Thus, to calculate an NMR spectrum, one^ constructs all of the possible 
basic product functions and sorts them into groups, each member of a 
group containing the same number of a and B spins. The energy of .sach 
of them is determined and used as the diagnonal element in a square matrix. 
The of f "diagonal elements of this matrix (H matrix) are energies of inter- 
action botv/een the members of the group. This matrix is diagonalized and 
these elements then contain the energies of the final spin functions. 
These final spin eigenf unctions are the columns of the matrix (U matrix) 
required to diagonalize the H matrix. Each element in a spin function 
represents the contribution of a basic product function to that spin 
function . 



Before continuing, it is necessary to define the value of a spin func- 
tion, as shown in equation 3, where s^ is as previously defined. 

3 F = P s. 

^ i=l ^ 

Nov/ the peaks observed in NMR spectra are transitions from one spin func- 
tion to anothor. Those transition Gnorgies and their intensities can be 
oalculrjted as outlined bolov;. Transitions are allowed between spin func- 
tions v/hose F^^ values differ only by 1. If this holds, the energy of the 
transition is the difference in energy of the final spin eigenf unctions , 
and the intensity of the transition is given by equation 4 where C^^ is the 
uth element in one spin function, is the vth element in the other spin 
function, and A-1 if the basic product functions represented by and 
differ by one spin. Otherwise, A==j3. 

4 I (Z I c C*A)^ 

U V U V 



8.0 EXAMPLE 

The following pages represent an actual printout. The spectra are repre- 
sentations of those displayed on the scope. The underlined data is that 
typed by the usr^r. 



LO NMRSIME,2 



NMRSIM 

COMMENTS : DJ-MO 1 

WAIIT PAPER TAPE I/0?Y 

WANT HIGH SPEED READER-PUNCH?y 

NO, OF SPINS= 5 

OFFSET & WIDTH 50^ 

CHEMICAL SHIFTS: 9j2f , 9j0 , 9j0 , 26^ , 260 

COUPLING CONSTANTS: ^ / Jg , 3 . 8 , 0 , 8 , 8 , 8 , 8 , JZI , 
BLOCK, U: 4 00, 3 



J 



L: 



MIN. INTENSITY :0 



DISPLAY LIST?Y 



INTENSITY 



ENERGY 
272.0.0 
264.00 
256.00 
248.00 
97.00 
90.00 
81.00 



27.00 
90.00 
9 9.00 
36.00 
104.00 
189.00 
86.00 



Pj 
0:" 
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OFFSET & WIDTH : J2f, 3jglj2f 

TAPE INPUT?Y_ 
BLOCK, U: 42jg,3, 

M0RE7N 
P: ^ 

NMRSIM 

COMMENTS: DEMO 1 CONT 

NO. OF SPINS= 2_ 

OFFSET & WIDTH ;^j_5g£ 

CHEMICAL SHIFTS: 265,280 

coupling; COiJSTAriTS: 5^ 

J:iLOCK., U: 4 40,3 



f~r-T-T r "T - T i--r -r- t"T"i — r T -y-r-r- i — r 'l i i i 

C^ 

COUPLING CONSTANTS: 1^ 

LLOCK, (J: 440_,3_ 
O: 

TAPE ir]PUT?N 
BLOCK, U: 440,3 



KORE?N 



T -T t ■ r-T 



D: 



9.0 ASSEMBLING NMRSIM 



The NMRSIM program is supplied to the user in both source and binary on 
the tape. To generate a new binary file if the source program is modified, 
use the following procedure. 



1. Load a DIAL-MS tape on unit 0, Load unit 1 with either 
a DIAL*-V2 or DIAL-MS system tape. If another tape unit 
is aval Labia, mount the tape containing the CLEARSYM and 
MMRSTM source programs there. IP only' two tape units 
are a va i. l.able^ p 1 ace the sourer; programs on unit 0 (with 
PIP if riocossary), in order to reduce the assembly time. 

2. Type ^ZK ^ to clear the binary V?orking Area of unit 1. 

3. Type ->AS CLLARSYM ^ . CLEARSYM is a two word program 
which produces a clean symbol table. 

^000 

SAVSYM 1 



4. Type >AS NMR2,0^ , Krror messages generated at this 
time should be ignored. Press the RETURN key to ter- 
minate the assembly after the errors have all been prin- 
ted to suppress printing of the symbol table. 



Typo >-7>eJ , This clears the binary Working AmsL on unit 1, 



Typo .LI IJMR3,0j. If no listing is desired, use the AS 
command. Any errors generated now are real and must be 
corrected . 



7, Type ->-LI NMR2,0^ . Because the symbol table produced is 
the same as the one generated in step 6, printing may be 
suppressed with the RETURN key after it has started to be 
printed. 



8. Type -SB NMR23,0} , This saves the binary output from the 
two previous assemblies 



9. Type -^ZeJ . 



10. Two versions of the Floating Point Packacjo are aipplied . 
If the raachine has the KAE option, NMRIL should bo used 
in the following steps and NMI^STMK v/ill bo generated. 
If the machine does nothave this option, NMRl should be 
used to generate NMRSIM. 
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11. 



Type -AB NMRl J . 



12, Type -AB NMR2 3,jef J . 

13. Type -^B NMRSIM.^PJ . A binary file of NMRSIM(E) has now 
been gencreitecl and tiB command ->L0 NMRSIM,0 ^. will cause 
load and. execute. 

10.0 INTERNAL DESCRIPTION 

The complete program listings are contained at the back of the manual. The 
following is a brief outline of the program; capital letters below refer to 

variables . 

1. Begins with text mode for input of commentary. 

2. Enter the number of spins and store in N. 

3. Set "switch" for punched output. 

4. Set up basic produce functions (BPFs). 

A, Calculate HARRY array v;hich contains relative addresses 
of the BPFs of different F value. 

z 



B. Calculate NUSE array. This array actually contains the 
BPFs. They are stored in groups, each group having an 
value one greater than the preceding group. 

5. Get sweep offset and sweep width. 

6. Get chemical shifts. 

7. Get coupling constants. 

8. Initialize SPEC array. 

A. This array eventually contains the displayed spectrum. 

a. It is initialized to jd except for calibration points 
where it is set to 4000^. 

o 

9. Calculate first H and U matrices. 

A. First H matrix is always a 1 X 1 which does not require 
di agona 1 i z a tion . 

B. First U matrix is always a 1 X 1 with its only element 
equal to 1. 

Set FLAG which is the negatlvo of the number of times the 
loop from 11 to 17 must be executed. 



10. 
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11. Set N2, calculate EN and UOLD arrays. 

A. Copy the trace of the H matrix into the EN array, 

B. Copy the U matrix into the UOLD array. 

C. Set N2 equal to the size of the H matrix just copied. 

12. Calculate the size of the next H matrix and put in into Nl. 

13. Determine which transitions will be allowed between the 
EPFs represented by UOLD and those that will be represented 
by U. This is stored in the array called TABLE. 

14. Initial izo II and U matrixes. 

A. Diagonal H elements are calculated by equation 1. 

B. Off-diagonal II elements are calculated by equation 2. 

C. The U a-nahrlx is initialized to all zeros except for 

ones on the diagonal. 

15. Diagonalize the 11 matrix and calculate the U matrix by 
the Jacobi method. 

16. Calculate transitions, 

A. Calculate transition energies by forming the differ- 
ences between the elements of the trace of the H mat- 
rix and the elements in the EN array. 

B. Calculate transition intensities from equation 4, 

C. Store appropriate values in SPEC array. 

17. Test FLAG to see if the loop is done. 

A. If loop is not done, go to 11, 

B. If loop is done, display spectrum, 

18. If ALT MODE is typed, interpret mnemonic and execute 
option. 
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NMRSIM 

i)EC»12-l^SA«D 



Dii[i\ii\ f'.qufpnicni C\)rpora{ ion maintains a continuous elToii to improve the quality and usefulnesj; of its 
pub]ii..,it ions. To do this effectively we need user feedback your critical evaluation of this manual. 

Ple.ise tjjnunent on tins nninuafs compIeter)ess, accuracy, orr^anj/ation. usability, and readability. 



D\ii vou find errors in this manual? , 



How can this nianual be improved'.'. 



^^^[^- ^'^^'"-^^-^ '<^-^P cu;>Ujmers informed of curreni DEC software and publications. Thus, tlie toliowmg pcriod- 
!c::jilv iUstnbuicd pubhcat ions are available upon request. Please check the appropriate boxes for a current issue of the 
p(ih!!ca!ion(s) desired. 

^ Manua l tJpdate. a quarterly collection of revisions to current software manuals. 

^ ^ B t jok'du-if, a bibliography of current software manuals. 

□ Prn|.'r,irn (.ihtary Pmcc \ a:U a iist of currently available software piogranis and manuals. 



i*k.a'.r describe your p<isi(ion. 



^''f^'^' Organization 

^^y^-^-^ ...^ Department 



( Stale „ Zip or Country. 



Do Not Tear - Fold Here and Staple 




HOW TO OBTAIN SOFTWARE INFORMATION 



AinHnincunscnls lor fk'w uml revised sollwurc, as wcH as pfograrnnniig notes, soStwarc in-obkMiis, anti tiocunicnla- 
iiun i..<trR-LUons arc pubhshod by Soflware [nioiniation Service in the lollowing ncwsjcflers. 



Digita! Soil ware News for the P[)f'-8 l-amiiy 
DiHital So! I ware News Tor the FDP-«i/15 Faniilv 
Pl)F-f)/fM)P- 10 Soil ware Builctin 



I hese tiewsjeUers cofjiaifi iiilornial loi! appHeable !<) s{>rtware available i'roni l>igitars Program Library. 



iMease <. otuplele Hte card below to {liaee yotir name on the riewsleUer inailing list. 



(,)ucstionN o! probieuis cutRerning Software slionkl be rcijoried to \\h: Sollware Speciaiist at your nearest DFC 
rejiional or district sales (jflice. In cases where no Software Specidlist is available, please send a Software Trouble 
Kejjort form with details ot the problem to: 

Software f iiforiDalion Service 
Dig! tall i\ u t p rn e n t ( 'orpor a 1 1 < ) n 
\4h Main Street, Bldg. 3-5 
Maynani, Massaciiusetts 0 ! 754 

[ hcsi; iornis. whici) are avaslable withorit cluirge from tlic Program l.iiirary, shotild be fully filled out and accompa- 
nied by teletype oulput as weil as lisln^gs or tapes of the user prograsii to facilitate a complete investigation. An 
answer will be sent to the itidsvidual and appropriate topics of general interest will be printed in the Jiewsletter. 

New and revised software and manuals, Software Trouble Report (ortns, and cumulative Software Manual Updates 
are avaiiable from the Program Library. When ordering, include the document number and a brief description of 
the program or manual rc(|uested. Revisions of programs and documents wit! be announced in the newsletters and 
a pnce list will be included twice yearly. Direct all inquiries and retjuests to: 

Program Library 
Digital h{|uM'^i^^^'^'s^ Corpora! i(.)n 
146 Main StreeL Bldg. 3^5 
.Maynard., Massacliusetts 0 1 754 

Inuii^u l:qu;prnent Computer Users Society (Dr.CUS) maintains a user Library and publishes a catalog of programs 
as well as the DHCUSCOPF magazine for its members and non-members who recjuest it. For further information 

please write to: 

DBCUS 

Digital i-qmpment Corporation 
146 Mam Street 

Maynard, Massachusetts 01754 



Send Digital's software newsletters to; 

Name 

Company Name 
Address 



My computcf is a PDP-8/1 □ 

LINC'8 D 
PDP-9 □ 
PDF- 10 □ 

My system .serial number is 



PDF-8/L □ 
PDP-I2 □ 
PDF- 15 D 

OTHER □ 

(if known) 



(zip code) 

Please specify 



Do Not Tear - Fold Here and Staple 
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Postage will be paid by: 



mill 



Digital Equipment Corporation 
Software Information Services 
146 Main Street, Bkig. 3-5 
Maynard, Massachusetts 01754 
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